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ABSTRACT 

Stellar models which incorporate simple diffusion or shear induced mixing 
are used to describe canonical extra mixing in low mass red giants of low and 
solar metallicity. These models are able to simultaneously explain the observed 
Li and CN abundance changes along upper red giant branch (RGB) in field low- 
met allicity stars and match photometry, rotation and 12 C/ 13 C ratios for stars in 
the old open cluster M67. The shear mixing model requires that main sequence 
(MS) progenitors of upper RGB stars possessed rapidly rotating radiative cores 
and that specific angular momentum was conserved in each of their mass shells 
during their evolution. We surmise that solar-type stars will not experience 
canonical extra mixing on the RGB because their more efficient MS spin-down 
resulted in solid-body rotation, as revealed by helioseismological data for the Sun. 
Thus, RGB stars in the old, high metallicity cluster NGC6791 should show no 
evidence for mixing in their 12 C/ 13 C ratios. 

We develop the idea that canonical extra mixing in a giant component of 
a binary system may be switched to its enhanced mode with much faster and 
somewhat deeper mixing as a result of the giant's tidal spin-up. This scenario 
can explain photometric and composition peculiarities of RS CVn binaries. The 
tidally enforced enhanced extra mixing might contribute to the star-to-star abun- 
dance variations of O, Na and Al in globular clusters. This idea may be tested 
with observations of 12 C/ 13 C ratios and CN abundances in RS CVn binaries. 

Subject headings: stars: chemically peculiar — stars: evolution — stars: interiors 
— stars: late-type — stars: rotation — globular clusters: general 
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1. Introduction 

Standard stellar evolution models assume that mixing only occurs in convective regions 
of stars. However, there is convincing evidence that a majority of low-mass (0.8 < M/M Q < 
2) stars are experiencing nonconvective extra mixing on the upper red giant branch (RGB) 
(e.g., Sweigart & Mengel 1979; Langer et al. 1986; Sneden et al. 1986; VandenBerg & Smith 
1988; Smith & Tout 1992; Charbonnel 1994, 1995; Charbonnel, Brown & Wallerstein 1998; 
Charbonnel & do Nascimento 1998; Gratton et al. 2000; Bellman et al. 2001; Grundahl et 
al. 2002; Denissenkov & VandenBerg 2003a; Smith & Martell 2003; Shetrone 2003). These 
red giants have a helium electron-degenerate core, a hydrogen burning shell (HBS) atop the 
core, and a thin radiative zone between the HBS and the bottom of the convective envelope 
(BCE) (Fig. 1). The most convincing observational evidence of mixing in red giants is the 
evolution of their surface carbon isotopic ratio 12 C/ 13 C and abundances loge( 7 Li), [C/Fe] and 
[N/Fe] 1 with increasing luminosity L along the RGB (circles in Fig. 2). The first abundance 
changes that occur on the subgiant and lower red giant branch (0.6 < log L/L Q < 1.0) are 
produced by the deepening BCE during the standard first dredge-up episode (Iben 1967). 
Here, the convective envelope engulfs the stellar layers in which the CN-cycle had altered 
the relative abundances of 12 C, 13 C and N while on the main sequence (MS). Convection also 
dilutes the envelope abundance of Li that has survived only in a thin subsurface layer where 
the temperature was less than ~2.5 x 10 6 K. At the end of the first dredge-up, the BCE 
stops to deepen and begins to retreat. It leaves a discontinuity of the chemical composition 
at the level of its deepest penetration. The evolutionary changes of the surface element 
abundances resume precisely at the moment when the HBS, advancing in mass inside red 
giants, erases the composition discontinuity left behind by the BCE. The commonly accepted 
interpretation of this is based on the following assumptions: (1) some extra mixing is present 
in the red giants' radiative zones; (2) it connects the BCE with the layers adjacent to the HBS 
(Fig. 1), where the CN-cycle is at work and Li is strongly depleted, thus dredging up stellar 
material enriched in N and deficient in C and Li; (3) at lower luminosities, extra mixing is 
shielded from reaching the vicinity of the HBS by a large gradient of the mean molecular 
weight (a /i-gradient barrier) built up by the composition discontinuity (Charbonnel, Brown 
& Wallerstein 1998; Denissenkov & VandenBerg 2003a). 

When the HBS meets and erases the /x-gradient barrier, the evolution of red giants slows 
down for a while. This results in prominent bumps in the differential luminosity functions 
of globular clusters, which have really been observed (e.g., Zoccali et al. 1999; Monaco et al. 



1 We use the standard spectroscopic notations: [A/B] = log (N(A)/N(B)) - log (N(A)/N(B))q, where 
N(A) and N(B) are number densities of the nuclides A and B; loge( 7 Li) = log(A( 7 Li)/A(H)) + 12. 
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2002; Cassisi et al. 2002; Salaris et al. 2002; Riello et al. 2003). Therefore, the luminosity 
at which extra mixing starts to manifest itself (at log L/L@ ~ 1.8-2.0 in Fig. 2) is called 
"the bump luminosity". It divides the RGB into its lower and upper part. In stars with 
M > 2M , the HBS fails to cross the /z-gradient barrier before the core helium ignition is 
triggered. In concordance with this, no manifestations of extra mixing have been reported 
in intermediate-mass red giants (Gilroy 1989). 

In the absence of a reliable theory of extra mixing in upper RGB stars, semi-empirical 
models are worth applying. In particular, a simple diffusion model has shown that in the 
majority of upper RGB stars the depth and rate of extra mixing do not seem to vary 
greatly from star to star. According to Denissenkov & VandenBerg (2003a), these quantities 
can be parameterized by any pair of correlated values within the close limits specified by 
A log T ps 0.19 and D mix ps 4 x K^cn^s -1 , to A log T ps 0.22 and D mix «8x 10 s cm 2 s _1 . 
Here, A log T is a difference between the logarithms of temperature at the base of the HBS 
and at a specified maximum depth of extra mixing M mix (Fig. 1), and D mix is a constant 
diffusion coefficient. 

Given that the majority of metal-poor upper RGB stars, both in the field and in globular 
clusters, experience extra mixing with almost the same depth and rate, Denissenkov & 
VandenBerg (2003a) have proposed to call this universal process "canonical extra mixing" . 
Although its physical mechanism is not identified yet, the phenomenon of Li-rich K-giants 
and its interpretation are likely to support a rotational mechanism. Indeed, most of the 
Li-rich giants are located above the bump luminosity (Charbonnel & Balachandran 2000). 
Their proportion among rapid rotators (v sin % > 8 kms -1 ) is ~50%. Compare this to ~2% 
of Li-rich stars among the much more common slowly rotating (vsm i < 1 kms -1 ) K-giants 
(Drake et al. 2002). Because for extra mixing driven by kinetic energy of rotation D mix oc v 2 , 
a tenfold increase of v would enhance D mix from its canonical value of ~ 10 9 cm 2 s _1 to 
~ 10 11 cm 2 s _1 . Denissenkov & Herwig (2004) have shown that precisely such enhancement of 
-Dmix is required to produce Li in K-giants via the 7 Be-mechanism (Cameron & Fowler 1971). 
Therefore, following Denissenkov & VandenBerg (2003a), we hypothesize that sometimes, 
when an upper RGB star gets spun up by an external source of angular momentum, the 
presumably rotational extra mixing in its radiative zone may switch from its canonical mode 
to an enhanced mode with much faster and somewhat deeper mixing. 

Denissenkov & Weiss (2000) have proposed that single Li-rich giants might have been 
spun up by engulfed massive planets. Alternatively, a red giant can be spun up by a close 
stellar companion. In this case, its rotation is accelerated by a drag force originating from 
viscous dissipation of kinetic energy of a tidal hump raised in its convective envelope by a 
companion star. As a result of this interaction, spin and orbital rotation of the red giant 
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get synchronized, and its initially elliptical orbit can become circular. Fig. 3 illustrates tidal 
synchronization in real stars. Here, we have plotted observational data of De Medeiros et 
al. (2002) on projected rotational velocities v sin? and orbital periods P of G and K giant 
components of field binaries with available orbital parameters. The values of v sini were 
multiplied by the factor 4/V « 1.27 that takes into account a random orientation of the 
stars' axes of rotation. Filled circles are binary systems with a circular or nearly circular 
orbit (those with eccentricities e < 0.10). Systems with eccentric orbits (e > 0.10) are 
represented by open circles. A set of theoretical curves is constructed by us under the 
assumption that Q = u, where Q and uj are the spin and orbital angular velocity of a red 
giant, for 4 values of the red giant's radius: R/R Q =5, 10, 20 and 40 (solid curves from 
left to right). The tidal circularization time is longer than the synchronization time (e.g., 
Hurley et al. 2002), therefore the giants represented by filled circles are most likely to have 
synchronized their rotation. 

The goal of this paper is twofold. First, in §3 we demonstrate that canonical extra 
mixing in solar-metallicity red giants in the open cluster M67 can be modeled using essentially 
the same parameters as those adjusted for metal-poor stars. After that, we attempt to model 
enhanced extra mixing in red giant components of tidally locked (corotating) binaries that 
are spun up as a result of tidal synchronization between their spin and orbital rotation (§4). 
In order to get a self-consistent picture of switching from canonical to enhanced extra mixing, 
in § 4 we make assumptions and use f2-profiles and parameters of extra mixing compatible 
with those discussed in §3. Finally, §5 contains a summary of our main conclusions along 
with a proposal of observational tests that could confirm or reject our models. 



2. Stellar Evolution Code with Rotation and Extra Mixing 

Evolution of rotating stellar models with extra mixing on the upper RGB is calculated 
from the zero-age MS (ZAMS) to the RGB tip with an upgraded version of the computer 
code used by Denissenkov & VandenBerg (2003a,b). The most recent update is the adoption 
of Alan Irwin's improved EOS 2 . In addition, the energy losses due to neutrino emission are 
now calculated with the code distributed by Itoh et al. (1996). Our stellar evolution code 
uses OPAL opacities (Rogers & Iglesias 1992) for temperatures above ~ 10 4 K, complemented 
by Alexander & Ferguson (1994) data for lower temperatures: both datasets assume solar 
abundances of a-elements. With the new EOS, we use a value of l com /H P = 1.75 for the 



2 We use the EOS code that is made publicly available at ftp://astroftp.phys.uvic.ca/pub/ 
irwin/eos/code/eos_demo_f ortran.tar .gz under the GNU General Public License. 
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ratio of the convective mixing length to the pressure scale height derived from a calibrated 
solar model. 

In our MS models with M > 1 M , some overshooting has been introduced beyond 
the formal convective core border determined by the Schwarzschild criterion by assuming 
that the ratio of the radiative and adiabatic temperature gradients (logarithmic and with 
respect to pressure) at the core border V ra d/V a d = 0.9 instead of 1. Using this approach, we 
succeeded in fitting our solar-metallicity 4x 10 9 yr isochrone 3 (Fig. 4, solid curve) for the high 
precision color-magnitude diagram (CMD) of M67 published by Sandquist (2004) (Fig. 4, 
open circles). We ignored any microscopic mixing, such as atomic diffusion, gravitational 
settling and radiative levitation. 

The treatment of extra mixing in the radiative zones separating the BCE from the HBS 
is the same as that described by Denissenkov & VandenBerg (2003a, see their equation (1)). 
Effects of shellular rotation on the internal structure and evolution of stars are treated in 
the same way as by Denissenkov & VandenBerg (2003b). The code employed in this work 
only considers the simplest cases of the angular momentum profile evolution: either (a) 
supporting the uniform (solid-body) rotation, Q(t,M r ) = £l(t), or (b) keeping the constant 
specific angular momentum in each mass shell, including convective regions, j(t,M r ) = 
|r 2 f2 = j(0,M r ). These two opposite cases correspond to a maximum possible and a zero 
efficiency of the angular momentum redistribution throughout a star. 

3. Canonical Extra Mixing in Single Red Giants of Solar Metallicity 

In this section, we present results of our new stellar evolution computations that, when 
being compared with recent observational data on photometry, rotation and chemical com- 
position of Population I (in the open cluster M67) and Population II (both in the field and in 
globular clusters) low-mass stars, allow to constrain the strength of rot at ion- induced shear 
mixing. 

3.1. Basic Equations and Assumptions in the Shear Mixing Model 

When studying the evolution of metal-poor stars with M 0.9 M q , Denissenkov & 
Tout (2000) and Denissenkov & VandenBerg (2003a) have proposed that canonical extra 



3 Necessary bolomctric corrections and T e g- color transformations were calculated with a code kindly 
presented to one of us (PAD) by Don VandenBerg. 
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mixing in upper RGB stars can be identified with turbulent diffusion induced by secular 
shear instability in their differentially rotating radiative zones and that it can be modeled 
using the coefficient of vertical turbulent diffusion derived by Maeder & Meynet (1996) 
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is the thermal diffusivity with k and Cp representing the opacity and the specific heat at con- 
stant pressure, respectively. The quantities 5 = — (dlnp/dlnT) p and (p = (d\np/d\np) PT 
are determined by the equation of state. 



Equation (1) gives simultaneously the extra mixing rate 
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and depth M m i x , the latter being determined as the coordinate of the first mass shell beneath 
the BCE where D v = 0, or 
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Although the depth given by equation (4) for reasonable fi-profiles is surprisingly close to 
that specified by the parameter A log T = 0.19-0.22 in our simple diffusion model, the rate 
calculated with equation (3) when / v = 1 is found to be too slow for canonical extra mixing 
(Denissenkov & VandenBerg 2003a). Therefore, we have to insert an enhancement factor 
f v 3> 1 into equations (1) and (3). 

For the low values of vsini usually measured in MS stars with M < 1M , rotational 
mixing, such as meridional circulation or turbulent diffusion, can reproduce the evolutionary 
changes of [C/Fe] in their upper RGB descendants only assuming that the specific angular 
momentum is conserved in each mass shell of these stars during their entire evolution from 
the ZAMS to the RGB tip (Sweigart & Mengel 1979; Smith & Tout 1992; Denissenkov & 
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Tout 2000; Denissenkov & VandenBerg 2003a). Even having assumed that, some differen- 
tial rotation in the initial models is still required. Following Denissenkov & VandenBerg 
(2003a,b), we set up such differential rotation in our ZAMS models assuming that the ratio 
of the equatorial centrifugal acceleration to the gravity is constant along the radius and that 
it is a small fraction of the critical ratio for the Roche potential used in our rotating stellar 
models; i.e., 

e{M r ) = = f £ e CTit = constant, (5) 

where e cr i t ~ 0.24. The ratio of the angular velocity near the center to Q at the surface 
obtained in this way is less than 10 for all of the ZAMS models considered below. Thus, we 
admit that the solid-body rotation of the inner Sun as revealed by helioseismological data 
(e.g., Chaplin et al. 1999) may not necessarily be applicable to all low- mass MS stars. The 
free parameter f £ is chosen so as to get theoretical values of v TOt close to projected rotational 
velocities of MS and subgiant stars. After that we evolve our models from the ZAMS to the 
bump luminosities keeping j constant in every mass shell. No extra mixing is allowed in the 
models until they will reach the bump luminosities. For the sake of simplicity, we ignore 
mass loss (in fact, it becomes important only in the very vicinity of the RGB tip). 



3.2. A Summary of Results for the Single Star Case 

Equations (4) and (3) with values of f2 and (<91nf2/<91nr) taken from the radiative 
zones of our bump luminosity models give approximately the same depths of extra mixing 
and diffusion coefficients that are only ~ 15-25 times as small as those constrained by our 
semi-empirical diffusion model. This conclusion is true not only for metal-poor low-mass 
stars, as was anticipated by Denissenkov & VandenBerg (2003a), but also for their solar- 
metallicity counterparts. To demonstrate this, we have computed the evolution of three 
rotating stellar models: the first two for the same mass M = 0.85 M and helium abundance 
Y = 0.24 but for two different heavy-element mass fractions Z = 0.002 and Z = 0.0005, 
and the third model for M = 1.35 M , hydrogen abundance X = 0.70 and the solar value 
of Z = 0.0188. 

The initial rotation profiles of the first two models have been specified by the parameter 
f £ = 0.0003. This results in surface equatorial rotational velocities v rot monotonously de- 
creasing from ~ 7kms~ 1 to ~4kms~ 1 during the models' MS lives. When these metal-poor 
models leave the MS and become subgiants, their rotation slows down to v rot < 4kms~ 1 due 
to a surface radius increase (Fig. 5, solid curve). Similar values of "true" mean rotational 
velocities have been measured in the MS turnoff and subgiant stars in the globular clusters 
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NGC 104, NGC 6397 and NGC 6752 by Lucatello & Gratton (2003) (squares with errorbars 
in Fig. 5). The third model evolves along the RGB of the 4 x 10 9 yr isochrone that matches 
the CMD of the old solar- metallicity open cluster M67 (Fig. 4, dashed curve). For this model, 
the parameter f £ = 0.00075 has been chosen so as the computed behavior of v rot along the 
rotational isochrone of the same age constructed using stellar models with M > 1.1 M and 
the same parameter f £ (Fig. 6, solid curve) to be in a qualitative agreement with the ob- 
served decrease of (4/ir)v sin i in the M67 stars whose evolutionary status changes from the 
MS to the subgiant branch (Fig. 6, circles). Unfortunately, at present we cannot discriminate 
between the cases of differential and solid-body rotation of M67 stars (compare solid and 
short-dashed curves in Fig. 6). However, the observational data in Fig. 6 seem to support 
our assumption of slow surface rotation (t> rot ~ 10 kms -1 ) of M67 stars on the early MS 
(compare solid and short-dashed curves with dot-long- dashed and long-dashed curves). The 
stars from Fig. 6 are also plotted in Fig. 4 as filled circles. Comparison of these two figures 
shows that the stars with (B — V)o ~ 0.55-0.60 in Fig. 6 have M < 1.35 M Q . Their masses 
are probably close to the M67 MS turn-off mass, which is M pa 1.25 M & . Only the objects 
with (B — V) > 0.9 in Fig. 6 are subgiants and lower RGB stars having M pa 1.35 M Q . 

In Fig. 7, vertical dotted and dashed line segments in panels d indicate depths specified 
by A log T = 0.19 and 0.22 in the third (M = 1.35 M , Z = 0.0188) unmixed model that 
is located immediately above its bump luminosity. Solid line segment points to the mixing 
depth calculated with equation (4) using an fi-profile extracted from the same stellar model. 
Note the proximity of the broken and solid line segments. 

Starting from the bump luminosity models, we continue our stellar evolution computa- 
tions. Rotational shear mixing described by equation (1) is now taken into account. To be 
more exact, we introduce a diffusion coefficient D mix = D v for D v > while letting D mix = 
for D v < 0. A new f2-profile required to recalculate D v is taken from the evolving rotating 
stellar models in each time step. As in the models below the bump luminosities, the specific 
angular momentum is still assumed to remain constant in each mass shell, including con- 
vective envelopes. Thus, our new extra mixing computations include for the first time both 
a plausible physical description for the diffusion coefficient D m[x , with its profile allowed to 
change along the whole length of upper RGB 4 , and a feedback effect of rotation and chemical 
mixing on the structure and evolution of red giants. In order to simulate the evolutionary 
variations of the surface 12 C/ 13 C ratio and abundances of Li, C and N in the field metal-poor 
upper RGB stars we had to increase by a factor of f v = 15-25 the value of D v given by 
Maeder & Meynet (1996) in our 0.85 M rotating models (Fig. 2). At the same time, the 



4 Denissenkov & VandenBerg (2003a) used a fixed D m i x -profile proportional to D v in the bump luminosity 
model. 
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12 C/ 13 C ratios in a couple of upper RGB stars and in a few clump stars in the cluster M67 
were reproduced by our 1.35 M & model only after D v had been multiplied by / v = 15 (Fig. 8, 
dashed curve). For comparison, dotted curve in Fig. 8 has been computed with the diffusion 
model using A log T = 0.22 and -D m i x = 8 x 10 8 cm 2 s _1 . Also, compare our Fig. 2 with 
Fig. 4 from Denissenkov & VandenBerg (2003a). The ratio between the enhancement fac- 
tors 25 and 15, adjusted for the low-mass red giant models of different metallicities, is even 
smaller than the one between our empirically constrained limits for the diffusion coefficient: 
Anix = 8 x 10 s cm 2 s" 1 and D mix = 4 x 10 8 cm 2 s _1 . 

Fig. 7 shows element abundance profiles, /i-gradients and diffusivities for two rotating 
stellar models with M = 1.35 M Q , Z = 0.0188 and f £ = 0.00075 located close to their bump 
luminosity. Compared to the case of metal-poor red giants, the only important qualitative 
distinction here is the additional hump in the /i-gradient profile at M r /M Q xs 0.2535 in the 
upper RGB model clearly seen in panel d. It is built up by H abundance variations that 
accompany the 12 C to 14 N transformation (see panel b). Because of this hump, the vertical 
turbulent diffusion cannot approach as close to the HBS as in the low-metallicity case. It 
can hardly reach a peak in the 13 C abundance profile (panel b). At the same location, the 
12 C depletion just sets in. This may explain why in the solar-metallicity case canonical extra 
mixing only produces modest changes of the surface 12 C/ 13 C ratio (Fig. 8). 

To summarize, multiplying the diffusion coefficient D v from Maeder & Meynet (1996) by 
a factor of f v = 20 ± 5, we can satisfactorily reproduce the observed evolutionary variations 
of the surface 12 C/ 13 C isotopic ratio and abundances of Li, C and N in the field metal-poor 
upper RGB stars (Fig. 2) as well as the scanty 12 C/ 13 C data for the upper RGB and clump 
stars in the old solar-metallicity open cluster M67 (Fig. 8). This can only be done under the 
assumption of differential rotation both in the low-mass MS ancestors of these evolved stars 
and in the convective envelopes of red giants. Otherwise, one has to search for an alternative 
mechanism of canonical extra mixing. 

Interestingly, our two models of canonical extra mixing (the diffusion model and the 
shear instability model) predict qualitatively different evolutionary changes of the surface 
Li abundance near the RGB tip. Whereas in the first model loge( 7 Li) remains very low, 
in the second model loge( 7 Li) resumes its post-first-dredge-up value and even exceeds it 
towards the RGB tip (compare upper panels in our Fig. 2 and in Fig. 4 from Denissenkov 
& VandenBerg (2003a)). This difference is due to a considerable increase of the diffusion 
coefficient D v oc K (equation (1)) in the second model caused by a growth of the thermal 
diffusivity K with the luminosity and the assumed constancy of D mix in the first model. The 
larger diffusion coefficient leads to some Li production via the 7 Be-mechanism. This new 
result can explain the higher Li abundances in clump stars compared to those in red giants 
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discussed by Pasquini et al. (2001). 

Now, we want to comment on the required enhancement of D v . We can formally comply 
with this requirement by choosing a value of the critical Richardson number Ri c larger than 
2 instead of taking / v >> 1. Indeed, the modern models of shear-driven turbulence do predict 
that Ri c should be at least 4 (Canuto 2002) or even 6.4 (Briiggen & Hillebrandt 2001) times 
as large as its classical value. If we accept this then the enhancement factor / v in equation 
(1) has only to be as large as ~ 20/6.4 m 3 for our model to reproduce the observations. It 
should also be noted that equation (1) takes into consideration only the largest turbulent 
eddies. Allowing for eddies from the whole turbulent spectrum would increase the coefficient 
D v in a similar way as the Full Spectrum of Turbulence model (Canuto & Mazzitelli 1991) 
generates much more vigorous convection than the Mixing Length Theory. 

As regarding the depth of shear mixing (4), it is directly constrained by the kinetic 
energy of differential rotation available to do work against the buoyancy force. It is our 
assumption of fl increasing with depth in the convective envelopes of low-mass red giants and 
in their MS ancestors that allows shear mixing to approach close enough to the HBS. In order 
to reconcile this assumption with the well established fact of the Sun's solid-body rotation, 
we refer to the recent work of Talon & Charbonnel (2004), according to which, in a solar- 
type star, angular momentum can be extracted from its radiative core by internal gravity 
waves generated in its convective envelope. The efficiency of this process should depend on 
the thickness of convective envelope. Therefore, we hypothesize that rotation-induced extra 
mixing may only work in those red giants whose MS ancestors had shallower convective 
envelopes than the Sun. Indeed, both our ZAMS models, the one with M = 1.35 M and 
Z = 0.0188 and that with M = 0.85 M and Z = 0.0005 possess much thinner convective 
envelopes than the solar ZAMS model (log M ce /M Q = —3.92, —2.32 and —1.62, respectively). 
Thus, we conjecture that metal-rich stars with M < 1.2 M Q might not experience canonical 
extra mixing on the RGB at all because their MS ancestors had thick surface convective zones 
and, as a result, nearly solid-body internal rotation, like the Sun. To test this hypothesis, it is 
necessary to determine 12 C/ 13 C ratios in upper RGB stars of an open cluster older than M67. 
A good candidate for these purposes is the open cluster NGC6791. It has [Fe/H] = +0.4 
and an age of (8.0 ± 0.5) x 10 9 yr (Chaboyer et al. 1999). 

4. Enhanced Extra Mixing in Giant Components of Tidally Locked Binaries 

Given that enhanced extra mixing is deeper than canonical one, it may penetrate the 
layers in outskirts of the HBS in which not only 12 C but also 16 is depleted in the CNO 
cycle and where 23 Na is produced in the reaction 22 Ne(p,7) 23 Na (Denissenkov & VandenBerg 
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2003a). Therefore, enhanced extra mixing in the extinct upper RGB stars that were slightly 
more massive (0.9 < M/M & < 2) than the MS turn-off stars in the present-day globular 
clusters (M rs 0.8- 0.9 M Q ) might have been one of the yet unidentified primordial sources 
of the global O-Na anticorrelation in globular cluster stars (Denissenkov & Weiss 2004). 

For the binary star case, we consider two basic stellar models: a solar-metallicity one 
with M = 1.7 M , and a metal-poor one with M = 1.0 M and Z = 0.0005. The first model 
has the mass and metallicity typical for both Li-rich K-giants and a sample of chromospher- 
ically active late-type giants — primaries of the RS Canum Venaticorum (CVn) binaries 
whose relevance to the problem of enhanced extra mixing will be discussed later. The sec- 
ond model represents a low-mass star that had completed its life in the past in a moderately 
metal-poor globular cluster. Both models will be "placed" into close binary systems with 
less massive stellar companions. Our main goal is to find out to what noticeable photometric 
and composition anomalies the tidal spin-up of these models may lead on the RGB. 



4.1. Description of Tidal Interaction in a Binary System 

4.1.1. Basic Equations 

If the initial separation between low-mass stellar components of a binary system is too 
small then the primary star can fill its Roche lobe before having reached the RGB tip. In 
this case, the radius of the red giant primary has an upper limit constrained by the Roche 
lobe radius R < Rl = E(q)a, where a is the binary semi-major axis, q = m/M is its mass 
ratio (secondary over primary), and, according to Eggleton (1983), 

°' 49 

^ ~ 0.6 + g 2 / 3 \n(l + q- 1 / 3 Y 

Dashed curve in Fig. 3 delineates the maximum surface rotational velocity of a red giant 
that fills its Roche lobe in a binary system with M = 1.7 M , q = 0.5, and an orbital period 
P, the latter being related to a via Kepler's third law 

(iH^ /3 (0 /3 < i+ ^ « 6 » 

where P^ is the period in days. We see that all stars but one lie below this curve. The only 
outlier is the G5III star HD 21018. It has v sini = 22.7 kms -1 , and the lithium abundance 
log£( 7 Li) = 2.93 (Hartoog 1978). It is most likely to be an intermediate-mass (M > 2.5 M Q ) 
star crossing the Hertzsprung gap on its way from the MS to the RGB region. Its location 
on the CMD (asterisk in Fig. 9) confirms our hypothesis. Thus, it may still preserve its 
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initial high Li abundance, while the large value of v sin % may be a signature of its previous 
fast rotation characteristic of the MS B-type stars. 

In our work, the tidal evolution of orbital and rotational parameters of binary star 
systems is modeled using differential equations derived by Hut (1981). Hut has elaborated 
upon the weak friction model of Alexander (1973), which assumes that the tidal humps 
in a gravitationally distorted rotating star have their equilibrium shape, as described by 
Zahn (1977), but with a constant time lag between the moment when the tidally deformed 
equipotential surface has been formed in the absence of dissipation and the current binary 
star configuration. For the sake of simplicity, we only consider the tidal evolution of primary 
stars possessing convective envelopes in binary systems with the zero eccentricity. In this 
case, only the spin rotation of the primary and the binary semi-major axis are affected by 
the tidal interaction, their changes being controled by the following equations: 

n = 3 ^iUJ ( } ' ( ) 

"i = -6^(1 + ,) a (8) 

at if V a J u 



where 



fMR 2 \ 1/3 
tconv ~ I I = 0.4304 



(M/M & )(R/R Q f ] ^ 



is a characteristic time of viscous friction, which approximately equals the convective time 
scale (Zahn 1977) (friction experienced by the tidal humps is produced by the convective 
turbulent viscosity in the envelope), L is the primary star's luminosity, r g is its radius of 
gyration that determines the total angular momentum I = r^MR 2 , Q and u are its spin and 
orbital angular velocities. 

Equations (7-8) contain the apsidal motion constant of the primary 

k 3 ~^)] fl0 ) 

: - 2{2 + v [a{R)]Y (iUj 

The structural function i](a) can be obtained by integrating the differential equation of 
Radau (e.g. Claret & Willems 2002): 

a^- + Q^(r] + l)+r](r]-l) = Q (11) 
da p 

with the boundary condition 77(0) = 0. Here, p(a) is the primary's density profile, and p(a) 
is the average density inside the equipotential surface (isobar) whose rotationally distorted 
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shape 5 is described as 

r = a[l — e(a)P 2 (cos 6)} ps a, 
where P 2 is the Legendre polynomial of order 2, 9 is the colatitude, and 

_ f](q) 2 a 3 
8 ~ 3GM P 

is the ratio of the equatorial centrifugal acceleration to gravity divided by 3, with Mp 
being the mass inside the isobar P. The quantities r\ and e are related to each other via 
i] = (dine /din a). 

The quantity e has already been used to set up differential rotation in our ZAMS models 
(equation (5)). In our solar metallicity ZAMS model with M = 1.7 M , the initial f2-profile 
is specified by the same parameter f £ = 0.00075 as in our single star model with M = 1.35 M & 
in §3. It results in the surface rotational velocity v TOt ~ 10 kms -1 . In fact, MS stars with 
M > 1.4 Mq have extremely thin surface convective zones, therefore their rotation is not 
thought to experience magnetic braking (or other type of braking whose efficiency depends 
on the thickness of convective envelope). Hence, our 1.7 M model should have rather been 
assigned v rot ~ 100 kms -1 instead of 10 kms" 1 , and it should have rotated as the solid body. 
However, this inconsistency does not matter when we study the evolution of this model star 
as a primary in a close binary system because, after having been synchronized with its orbital 
rotation on the RGB, the primary will acquire a surface rotational velocity by far exceeding 
that it would have if it were a single star. So, the information about its previous surface ro- 
tation will be completely lost. What really matters is that our differentially rotating ZAMS 
model has approximately the same value of Q within its inner region, that will be occupied 
by the radiative zone on the upper RGB, that it would have if it rotated as a solid body with 
v TOt ~ 100 kms -1 . Indeed, in Fig. 10 long-dashed horizontal line shows the uniform rotation 
profile with t> rot ps 100 kms -1 , while dot-long-dashed curve presents the ^-profile obtained 
using equation (5) with f £ = 0.00075. Both profiles almost coincide at M r < 0.4 M Q . There- 
fore, under the assumption of the conservation of the specific angular momentum j in each 
mass shell, including convective regions, radiative zones of the single RGB models that have 
started their evolution on the ZAMS with these two different rotation profiles (the uniform 
one with v Iot 100 kms -1 , and the differential one with v Iot ps 10 kms -1 ) will possess simi- 
lar differential rotation. Dot- short- dashed curve in Fig. 10 demonstrates the rotation profile 
in the bump luminosity model obtained from the initial dot-long-dashed profile assuming 



5 According to Kippenhahn & Thomas (1970), within a high accuracy the structural changes of stars in 
tidally locked binaries are caused by their rotation and not by tidal forces. Therefore, we do not take into 
account the tidal effects on the primary's structure and evolution. 
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j(t, M r ) = j(0, M r ) for < M r < M. For comparison, short-dashed horizontal line gives the 
uniform fi-profile resulted from the initial long-dashed one assuming the solid-body rotation 
of the whole star during its entire evolution. In all the cases, the total angular momentum is 
conserved. Since the bulk of the red giant's angular momentum is contained in its convective 
envelope, the short-dashed line also approximates quite well the envelope rotation (on the 
right of vertical dotted line) in the alternative case, that we did not consider, in which j is 
only conserved in radiative regions, while Q is maintained constant in convective regions. In 
that case, as we mentioned, the radiative zone would rotate much slower than in the case 
represented by the dot-short-dashed curve (compare the values of f2 on the short-dashed and 
dot-short-dashed profiles at the BCE). 

Before "placing" our models into binary systems, we have computed their single stellar 
evolution. For the metal-poor model (M = 0.85 M , Z = 0.0005), the initial internal Q- 
profile has been specified by the same parameter f £ = 0.0003 as in §3. During the single 
stellar evolution computations, we were tabulating the quantities R, k and r g as functions 
of age and luminosity. On the RGB, the ratio k/r^ has been found to decrease very slightly, 
from about 0.5 to nearly 0.2 for both models. Hence, from equations (7) and (9) we can 
estimate the primary's synchronization time (the time after which Q m uo) as 

tp ( a \ 6 

^ 0.4304 

Furthermore, for values of q > 0.5, the product in front of the last term in equation (12) does 
not deviate much from unity during the entire evolution of our models from the subgiant 
branch to the RGB tip. Therefore, for binaries with the mass ratio 0.5 < q < 1, one can 
use the simpler estimate t syn ~ (a/R) 6 ji. Since the low-mass stars spend a time of the 
order of 10 8 yr on the RGB, only primary stars in the binaries with a < 20 R are expected 
to get synchronized on this evolutionary stage. If we want our primary to be tidally spun 
up already on the lower RGB, where R < lOi? , we should only consider the binaries with 
a < 200 Rq, or, in other words, only those with Pd < 300 (equation (6)). 

4-1.2. Choosing a Rotation Profile in the Convective Envelope 

Unfortunately, we do not know how rotation, convection and tidal friction interact with 
each other in the red giant's envelope. Even 3D-hydrodynamical simulations will hardly 
be able to model such complex situation with confidence in the nearest future. Therefore, 
we cannot do anything else but to make some plausible assumptions about the outcome of 
this interaction and to use as many as possible relevant observational data to constrain our 
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assumptions. We should also choose rotation profiles and parameters of extra mixing in our 
binary red giant models consistent with those employed in § 3 for the single star case. 

In Fig. 11, we have plotted fi-profiles in our solar-metallicity 1.7 M lower RGB models 
that started their evolution (on the ZAMS), either as a single or a binary star, with the same 
rotation parameter f e = 0.00075, and so with the same value of t> ro t, but for which different 
assumptions about the rotation law in their convective envelopes were made. All the models 
except that plotted with short-dashed curve have reached nearly the same evolutionary point 
immediately below the bump luminosity. Dot-dashed curve gives the f2-profile in the single 
star model. Long-dashed, solid and short-dashed curves show internal rotation of red giants 
in binaries with the same initial parameters: q = 0.5, a = 80 R®. Dotted curve has q = 0.5 
and a = 50 R Q . 

The long-dashed curve is obtained under the assumption that, after the corotation has 
been achieved (Q = uj) as a result of tidal interaction in the close binary, the whole convective 
envelope of the red giant rotates as the solid body, i.e. Q(r) = uo for -Rbce < r < R. However, 
in this case the radiative zone would rotate slower than in the single star case (compare the 
parts of the dot-dashed and long-dashed curves just beneath the BCE, at log R/R & < —0.2). 
Given that on the RGB hydrogen rich material is slowly moving from the convective envelope 
to the helium core through the radiative zone and the specific angular momentum is assumed 
to be conserved in this zone, we anticipate that, in the more evolved models, the f2-profile 
developed from that represented by the long-dashed curve will be parallel to the dot-dashed 
curve (like the short-dashed curve showing the advanced evolution of the solid curve) but 
it will entirely lie below it. In that case, since we assume that D mix oc Q 2 (e.g., equation 
(1)), the rate of extra mixing would be reduced compared to the canonical case, in spite 
of the red giant's faster surface rotation. This seems to be at odds with the reported very 
high percentage of Li-rich giants among rapidly rotating stars (Drake et al. 2002) and the 
necessity to enhance D mix from ~ 10 9 cm 2 s _1 to ~ 10 11 cm 2 s _1 in order to enable the 7 Be- 
mechanism (Denissenkov & Herwig 2004). Therefore, we reject the f2-profile represented by 
the long-dashed curve as implausible. 

In § 3, we have made the assumption that even convective envelopes of single stars rotate 
differentially, with the specific angular momentum conserved in each mass shell inside them, 
during the entire evolution from the ZAMS. Without that assumption rotational mixing 
driven by the secular shear instability in the radiative zones of upper RGB stars would be 
too weak to be identified with canonical extra mixing. At the same time, we admit that 
tidal interaction in a close binary system with a red giant component can enforce solid- 
body rotation of at least an outer part of its convective envelope. Making an appeal to 
observational data, we first notice that some of the chromospherically active late-type giants 
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in the RS CVn binaries (for instance, those shown as filled circles in Fig. 9) have almost 
equal orbital and rotational periods. This means that their surface layers have already been 
spun up by the tidal drag force to the orbital angular velocity, Q(R) = u. Although we 
know nothing about their internal rotation in the convective envelope below the surface, 
it is naturally to think that the same tidal force induces corotation of their subsurface 
convective layers as well, at least down to some depth R cor , i-e. Q(r) = u for R COT < r < R, 
where -Rbce < -Rcor < R- We have seen that the assumption of R COT = Rbce (long-dashed 
curve in Fig. 11) appears to be wrong if we want to explain canonical extra mixing and the 
phenomenon of Li-rich giants assuming that D mix oc Q 2 in both cases. 

The idea that R COT = R while j = constant in the convective envelope is not good either 
because in that case rotation of the red giant's radiative zone near the HBS £l(R c ) could 
easily become supercritical. Indeed, if in the whole radiative zone j = constant as well, then 



The hydrogen burning shell always stays at r m R c m 0.02 R e , and the helium core mass 
in stars located close to the bump luminosity is M r ks 0.3 M , which results in e(r)/e crit « 
0.12 (i?/i? ) 2 (f rot /lOkms -1 ) 2 . Thus, in a binary RGB star that has reached a corotation at 
the surface with t> rot = 10 kms _1 (Fig. 3 shows that this velocity is still less than the maximum 
possible value), rotation near the HBS could become critical already for R « 2.9 R Q , i.e. far 
below the bump luminosity. This disagrees with the presence of corotating binary red giants 
close to and even above the bump luminosity. 

Equations (7-8) were derived assuming for simplicity that Q = constant throughout 
the convective envelope. This may be a good assumption for thin convective envelopes but 
an RGB star has a very thick convective envelope with i?/i?BCE ^ 1- We can divide it 
into a number of thin spherical layers and consider them as rotating independently of one 
another. This approach is consistent with the assumption of differentially rotating convective 
envelopes made for single giants. If the tidal drag force had not varied along the radius inside 
the convective envelope then it would have naturally been to think that all the layers get 
synchronized at the same time. However, the transverse component of the tidal force strongly 
depends on r (Hut 1981), 



Therefore, the deeper convective layers are most likely to get synchronized much later than 
the layers near the surface. 

In our work, f2-profiles in convective envelopes of RGB model stars in close binary 
systems are computed as follows: (i) using the values of R, k and r g tabulated as functions 
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of age and luminosity during the computations of the single stellar evolution, we first solve 
equations (7-8) for a specified set of parameters q and a; this gives us dependences of v TOt 
on L (or My), and fl/u on R, like those plotted in Figs. 12 and 13 for our solar metallicity 
1.7 M Q model in a binary system with the mass ratio q = 0.5 (dashed curves are obtained for 
a = 50 Rq (P ~ 25 days) while solid curves are for a = 80 _R (P ~ 50 days); in Fig. 12, solid 
and dashed curves are computed (i.e. equations (7-8) are solved) starting with a subgiant 
model that has R = 3R Q , assuming that initially Q = 0, while dot-dashed and dotted curves 
- for f2 = 2.5 x 10 _5 rads~ 1 (this value of f2 extrapolates to v rot « 100 kms" 1 back on the 
ZAMS); (ii) from Fig. 13, we can read the star's radius R COT at which Vt/uj = 0.9; for the sake 
of simplicity, we assume that precisely at the moment when R = R cor the star's surface comes 
to corotation with the binary orbital motion; we find that R COT ~ 3.9 Rq and R cor m 6.7 Rq 
for a = 50 R® and a = 80 Rq, respectively (Fig. 13); (Hi) now, we "place" our model star 
into a binary system; until the star's radius R < R cor , its evolution is not considered to differ 
from the single star case, in particular, the whole star including its convective envelope is 
assumed to rotate differentially with j(t,M r ) = j(0,M r ); however, as soon as R > R cor , 
the solid body rotation f2(r) = uj of the outer part of the convective envelope, in which 
-Rcor < r < i?, is enforced. 

Our algorithm can create a faster rotating radiative zone in a red giant component of a 
wider binary in spite of the fact that a closer binary has a higher surface rotational velocity 
(compare solid and dotted rotation profiles in Fig. 11). Indeed, if in the case of a tidally 
spun-up RGB star the specific angular momentum is still conserved in each mass shell with 
R c < r < -Rcor, then the angular velocity in the radiative zone will depend on the radius as 
fi(r) = uo(R COT / r) 2 , where 



Using this equation, for our solar metallicity red giant model in the binaries with a = 50 R & 
and a = 80 R Q we find ^ = 0.2834 (R COI = 3.9 Rq) and uo 2 = 0.1400 {R cor = 6.7 Rq), 
respectively (both velocities are expressed in units of 10 _5 rads _1 and plotted in Fig. 10). 
Hence, at the same location in the radiative zone, the red giant in the binary with a = 80 Rq 
rotates 1.5 times as fast as the one in the binary with a = 5Oi? , in spite of the fact that 

UJ 2 <UJ\. 
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4.2. Possible Manifestations of Enhanced Extra Mixing 

4-2.1. RS CVn binaries 

As primary constraints on possible manifestations of enhanced extra mixing in giant 
components of tidally locked binaries, we have chosen the latest observational data on the 
CMD locations, rotational and orbital periods, rotational velocities, and surface element 
abundances for a sample of late-type red giants in the RS CVn binaries published by Fekel 
et al. (2002), Morel et al. (2004), and Fekel & Henry (2005). These red giants are chromo- 
spherically active stars whose activity is thought to be due to their tidal spinning up that 
assists in generating strong magnetic fields via convective dynamo. In order to be sure that 
all stars in the sample have achieved corotation, we select only those of them that have 
almost equal spin and orbital periods. Thus selected objects are shown as filled circles in 
Fig. 9. In the same figure, we have plotted evolutionary tracks for single solar-metallicity 
stars with masses 1.3, 1.5, 1.7, 2, and 3M (dashed curves). Looking at Fig. 9, one cannot 
help but noticing that all 8 of the selected RS CVn binaries are located on the lower RGBs. 
This striking photometric peculiarity has been emphasized by Morel et al. (2004). To extend 
the size of the sample of late-type red giants in tidally locked binaries, we have added a num- 
ber of K-giants in binaries with circular orbits from Fig. 3 to our selected RS CVn systems 
(open circles in Fig. 9). We have used the Hipparcos parallaxes and Tycho magnitudes from 
the online catalogue archive.ast.cam.ac.uk to estimate the absolute magnitudes My of 
the additional stars. The data plotted in Fig. 9 show that the majority (15 out of 17) of the 
observed low-mass (M < 2 M ) binary stars, that are expected to experience extra mixing 
on the RGB and that are most likely to have synchronized their spin and orbital rotation, 
are still located on the lower RGB. We will show that this photometric peculiarity of the RS 
CVn and circularized binaries can be explained by rotational effects in the radiative zones 
of their tidally spun-up red giant components. 

Solid curve in Fig. 9 is the evolutionary track of our solar-metallicity 1.7 M model 
placed into a binary system with q = 0.5 and a = 80 R & . As explained before, until the 
model's radius R stays below the corotation radius R cor = 6.7 R Q (for My < 1.68 on the 
RGB), its evolution is not considered to differ from that of the single star. In particular, the 
surface rotational velocity v rot decreases in a way similar to that depicted by the dot-short- 
dashed curve in Fig. 6. However, as soon as R exceeds -R CO r, the outer part of the model's 
convective envelope, at R cor < r < R, is enforced to rotate like the solid body with the 
respective orbital velocity f2(r) = uj 2 (see Fig. 10). It means that from the moment when 
R becomes equal R COI on the RGB (My 1.68) the surface rotational velocity v Iot of our 
model starts to increase following the thin solid curve in Fig. 12. After this moment, while 
the surface angular velocity of our binary model's single counterpart continues to decrease 
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due to the conservation of the total angular momentum and the red giant's expansion, the 
outer convective layers of our binary model keep rotating at the same angular velocity uo 2 
and its surface rotational velocity v Iot = Ruj 2 increases with R. By the bump luminosity, the 
ratio of uj 2 to Q(R) for the single model grows up to ~ 10 (compare the dot-dashed and solid 
curves at the maximum r in Fig. 11). At the BCE, this ratio is somewhat smaller because in 
the convective envelope of the single model Q(r) increases with decreasing r all the way from 
R to -Rbce, while in the binary model fl(r) is maintained equal to u 2 from the surface down 
to the depth R cor . Beneath the BCE, the same difference between the angular velocities 
for the binary and single star case that has been settled at the BCE is slowly propagating 
through the radiative zone, like a wave, during the subsequent stellar evolution (see the solid 
and short-dashed curves in Fig. 11). 

Usually, when the HBS crosses the chemical composition discontinuity left behind by 
the BCE at the end of the first dredge-up, a low-mass star makes a tiny zigzag on the CMD, 
like that outlined by a small square in Fig. 4. However, in our case of the solar-metallicity 
1.7 M Q binary star the zigzag is found to be unusually big (the solid curve between My pa 0.7 
and My pa 1.5 in Fig. 9). This is a new result which is entirely due to rotational effects. 
Indeed, in panel a in Fig. 14 the corotation radius is varied from 2 to 8R & , thus accelerating 
rotation of the radiative zone in the bump luminosity model of our 1.7 M & binary star, while 
panel b demonstrates how the zigzag's extent responds to the variations of R cor . These 
test computations have been done without extra mixing, just with rotational effects on the 
stellar structure and evolution included according to Denissenkov & VandenBerg (2003b). 
The region in the CMD in which our binary star makes the zigzag is comparable by its 
size with the area occupied by the majority of the observed tidally locked binaries (Fig. 9). 
Therefore, we surmise that it is the extended bump luminosity zigzags produced by tidal 
spinning up of the convective envelopes and underlying radiative zones that are responsible 
for the fact that the late-type giant components of the synchronized binaries predominantly 
reside on the lower RGB. This hypothesis receives a further support if we compare the 
evolutionary times or differential luminosity functions of the single and binary model stars 
in the vicinity of the bump luminosity. Whereas the time spent by the single star (thick 
curve in Fig. 15) between My = 2 and My = is 109 million years, it takes 146 million 
years for the binary star (thin curve in Fig. 15) to make the zigzag on the CMD when its 
magnitude My first increases from 0.7 to 1.5 and then decreases back to 0.7. Besides, the 
bump luminosity itself is shifted by AMy pa 0.8 towards the subgiant branch in the binary 
star case, which might explain the fact that the low-mass red giants in tidally locked binaries 
"prefer" to reside on the lower RGB. 

For 7 out of our 8 selected RS CVn binaries, measured values of v sin i are available 
in the cited papers. 4 binaries have orbital periods close to 50 days, for the rest 3 of them 
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P 7n 25 days. In most cases, an estimated mass of the primary component lies between 1.5 
and 1.7 M®, and its metallicity —0.09 < [Fe/H] < 0.12. That is why we have chosen the solar- 
metallicity 1.7 M model star placed in binaries with a = 50 R & (Pd ~ 25) and a = 80 P© 
(Pd ~ 50) as one of our basic binary models. Note that for Pd = 25 and Pd = 50, our 
binary 1.7 M© model with q = 0.5 will start to transfer matter to its companion by Roche- 
lobe overflow at My ~ —0.03 and My —0.65, respectively, i.e. well above the observed 
location of the RS CVn binaries in the CMD (Fig. 9). Therefore, their scarcity on upper 
RGBs can hardly be attributed to a dissipation of their envelopes on a short (dynamical) 
timescale in a common envelope event that follows the Roche-lobe overflow. Fig. 12 shows 
that, within observational errors, our solutions of equations (7-8), that control the tidal 
evolution of the red giant's angular velocity and the binary's semi-major axis, conform with 
the observed locations of the selected RS CVn binaries on the My versus i>sini plane. 

The tidal spin-up of the radiative zone may transform extra mixing from its canonical to 
an enhanced mode, with the mixing depth and rate increased proportionally to Q 2 according 
to equations (3-4). Dot-dashed curves in panels b and d in Fig. 16 demonstrate how the 
diffusion coefficient D v from equation (1) with f v — 1 changes when the envelope material 
spun up by the tidal force eventually arrives at the HBS in our 1.7 M Q upper RGB 6 binary 
model. Comparing panels c and d, we find that the tidally enforced enhanced extra mixing 
can dredge up some fresh Na that is being synthesized from 22 Ne in the HBS. This finding 
is not surprising because it has already been predicted for spun-up upper RGB stars by 
Denissenkov & VandenBerg (2003a). What is really surprising is that Na overabundances of 
the right magnitude have recently been reported in 5 out of our 8 selected RS CVn binaries 
by Morel et al. (2004). These are plotted in Fig. 17 against the Li abundances measured 
by the same authors. They have applied NLTE corrections to both [Na/Fe] and loge( 7 Li). 
They have also claimed that "chromospheric heating appears insufficient to account for the 
anomalously high Na abundances" . 

We have included the effect of tidal spin-up on D v in our stellar evolution computations. 
The resulting enhanced extra mixing produces dashed curve in Fig. 17. This curve corre- 
sponds to the stellar evolution that starts at the bump luminosity, includes the rotationally 
extended zigzag along the lower RGB and ends well above the bump luminosity. However, 
since most of the tidally synchronized binary red giants (including all 8 of our selected RS 
CVn binaries) in Fig. 9 are still located in the zigzag region, only the part of the dashed curve 
that corresponds to the evolution before the end of the zigzag should be compared with the 



6 We still call "the upper RGB models" those spun-up binary red giants in which the HBS has erased the 
composition discontinuity left behind by the BCE at the end of the first dredge-up in spite of the fact that 
they can reside on the lower RGB for quite a long time thereafter. 



- 21 - 



observational data points (a region between dotted line segments in Fig. 17). By the end of 
this part, [Na/Fe] reaches 0.36 and log£( 7 Li) accumulates only -0.6 after its initial drop by 
more than 3 orders of magnitude. These values are smaller than the observed abundances. 

A much better result (solid curve in Fig. 17) is obtained with the constant diffusion coef- 
ficient -D mix = 10 11 cm 2 s _1 (dotted curve in panel d in Fig. 16). In this case, we have actually 
made three simplifying assumptions: (i) that canonical extra mixing has operated with the 
constant diffusion coefficient D mix = 10 9 cm 2 s -1 (dotted curve in panel b in Fig. 16), which 
is close to the semi-empirical values estimated by Denissenkov & VandenBerg (2003a); (ii) 
the tidal synchronization of our binary red giant model has increased Q in its radiative zone 
exactly by the factor of 10, which is even slightly less than the ratio between the f2-profiles 
for the binary and single star case seen in Fig. 11; (Hi) like D mix from equation (3), our uni- 
form diffusion coefficient scales as D mix oc Q 2 . It is important to note that, notwithstanding 
these approximations, the mixing depth has yet been determined by equation (4) because we 
believe firmly that it is constrained by the kinetic energy available from differential rotation 
of the radiative zone. 

The difference between the dashed and solid curve in Fig. 17 is mainly due to the 
behaviour of the Li abundance, which can be explained by the same reason why loge( 7 Li) 
changes differently in a single star on the upper RGB towards the RGB tip in the secular 
shear instability and simple diffusion model of canonical extra mixing (compare upper panels 
in our Fig. 2 and in Fig. 4 from Denissenkov & VandenBerg (2003a)). In the first case, 
-Dmix is proportional to the thermal diffusivity K (equation (3)) that strongly decreases 
with the luminosity when our tidally spun-up red giant descends along the extended bump 
luminosity zigzag (solid curve in Fig. 9). Thus, on the one hand, the tidal spinning up tends 
to enhance -D mix because D mix oc Q 2 but, on the other hand, the effect of fast rotation on 
the red giant's internal structure causes its luminosity, and hence the value of K in the 
radiative zone, to drop, so the resulting mixing rate -D mix oc Vt 2 K turns out to be too slow, 
compared to the case of constant D mix = 10 11 cm 2 s _1 , to produce and support (via the 
7 Be-mechanism) as high Li abundance as in the RS CVn binaries (Fig. 17). Of course, 
if, in concordance with their apparent residence on the lower RGB (Fig. 9), the RS CVn 
binaries have not reached their bump luminosities yet then their surface Li abundances 
are in good agreement with those predicted by the standard theory for the first dredge- 
up (solid curve at [Na/Fe] < 0.05 in Fig. 17). However, in that case, RS CVn binaries 
would not have had as large Na overabundances as those reported by Morel et al. (2004). 
A spectroscopic determination of CN abundances in the giant components of the RS CVn 
binaries could verify our hypothesis of enhanced extra mixing in these stars. Indeed, if all the 
assumptions made are correct then, according to our computations, these stars should have 
12 C/ 13 C w 3.9, [C/Fe] ^ —1.8 and [N/Fe] ^ 0.63 instead of the standard post-first-dredge-up 
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values of 12 C/ 13 C^ 25, [C/Fe] « -0.17 and [N/Fe] « 0.30. 

4-2.2. Binary Red Giants in Globular Clusters 

The problem of Na overabundances in RS CVn binaries and its possible solution pro- 
posed by us may be related to the long-standing problem of star-to-star abundance variations 
of C, N, O, Na, Mg and Al in globular clusters. At present, there are no doubts that a size- 
able component of those variations originated from hydrostatic hydrogen burning in more 
massive stars that had completed their lives in the past. Indeed, firstly, there are correlations 
between overabundances of N, Na and Al and deficits of C, O and Mg as if all of them re- 
sulted from simultaneous operation of the CNO, NeNa and MgAl cycle. Secondly, the same 
abundance variations that are seen in red giants are also found in subgiant and MS stars, 
the latter evidently having too low central temperatures to produce these variations in situ. 
Hence, some of the globular-cluster MS stars must have accreted large amounts of material 
(or they have entirely been formed of it) processed in H-burning that had been ejected by 
now extinct stars into the interstellar medium of globular clusters. 

Figs. 18 and 19 give nice examples of the O-Na anticorrelation and Na-Al correlation 
for large samples of red giants in the globular clusters M3 (circles) and M13 (triangles) 
recently updated by Sneden et al. (2004). If we take [0/Fe] init = 0.3, [Na/Fe] init = -0.2 and 
[Al/Fe]i n it = —0.2 as the initial abundances in those M3 and M13 MS stars that were not 
polluted by their more massive cluster-mates, and [0/Fe] acc = —1.0, [Na/Fe] acc = 0.7 and 
[Al/Fe] acc = 1.2 as the most extreme abundances in the material accreted by some of the 
cluster MS stars then both correlations can be considered as simple mixtures of a (1 — at- 
traction of the material with the initial abundances and an x-fraction of the (accreted) 
material with the most extreme abundances (solid curves in Figs. 18 and 19; thick dots on 
the curves mark the values of x from to 1 with the increment 0.1). So, the next question 
is what H-burning stars were capable of producing those extreme abundances of O, Na and 
Al. There are two alternative answers to this question: (1) intermediate-mass (4-6M ) 
asymptotic giant branch (AGB) stars (e.g. Ventura & D'Antona 2005), or (2) upper RGB 
stars that had been more massive (say, 0.9 < M/M & < 2) than the present-day globular- 
cluster MS turn-off stars and that had experienced enhanced extra mixing caused by their 
tidal spin-up in binary systems (Denissenkov & VandenBerg 2003a; Denissenkov & Weiss 
2004). 

Since the first alternative is not without its shortcomings (e.g., when destroying O 
in intermediate-mass AGB stars the (hydrogen) hot-bottom burning depletes 24 Mg even 
stronger (Denissenkov & Herwig 2003); it also keeps [C/Fe] > —0.5 (Denissenkov & Weiss 
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2004); besides, the third dredge-up in AGB stars should increase the total abundance of 
C+N+O (Fenner et al. 2004); none of these theoretical predictions is supported by observa- 
tions), the second idea is worth developing. Pursuing this goal, we have computed the evo- 
lution of the 1 M Q star with the helium mass fraction Y = 0.24 and Z = 0.0005, which gives 
a metallicity [Fe/H] close to those of both M3 ([Fe/H] = -1.57) and M13 ([Fe/H] = -1.54). 
In order to simulate tidally enforced enhanced extra mixing, this model star has been placed 
into a binary system with q = 0.3 and a = 50 R & . We used the same computational method 
and approximations as in the case of RS CVn binaries. The initial rotation parameter was 
f £ = 0.0003, as in our 0.85 M Q single star model with Z = 0.0005, and the enhancement 
factor in equation (3) was / v = 20. The computed evolution of the O and Na surface abun- 
dances is shown in Fig. 18 with dashed curve. In order to produce enough Al (dashed curve in 
Fig. 19) we had to make the same additional assumptions as in Denissenkov & Tout (2000), 
namely: (1) the initial abundance of the 25 Mg isotope was increased to [ 25 Mg/Fe]i n i t = 1.2 
compared to the scaled solar ratio, and (2) the reaction 26 Al g (p,7) 27 Si was sped up by the 
factor of 10 3 compared to its rate from Angulo et al. (1999). 

Comparing the binary red giant evolution with that of a single red giant with the same 
mass and composition, we have found that the tidal spin-up has accelerated rotation of 
the radiative zone by a factor of ~ 20. Taking into account intrinsic imperfections of our 
description of the tidal spin-up, we have also considered a test case in which rotation of 
the red giant's radiative zone has artificially been maintained 25 times as fast as in the 
single red giant model. For this test case, the correlated O, Na and Al abundance variations 
are depicted with dotted curves in Figs. 18 and 19. In both cases, the abundances of O, 
Na and Al naturally evolve from their initial values of [0/Fe]i n i t = 0.3, [Na/Fe]i n i t = —0.2 
and [Al/Fe]i n it = —0.2 close to the most extreme values supposedly present in the material 
accreted by globular-cluster MS stars. 

It is important to note that environmental conditions in globular clusters are likely 
to be appropriate for our binary star scenario to contribute to the star-to-star abundance 
variations. Firstly, frequent single-binary and binary-binary stellar encounters work toward 
decreasing the semi-major axes of hard binaries (Hut et al. 1992). In particular, Beer & 
Davies (2004) have recently shown that after ~20 of such encounters the initial loga-flat 
binary population transforms into a gaussian-like distribution with a peak at a ~ 100 R Q . 
Secondly, binary red giants with a < 100 R Q will definitely fill their Roche lobes before 
reaching the RGB tip. Therefore they will probably lose their convective envelopes very 
quickly during a common-envelope phase (Beer & Davies 2004). According to our hypothesis, 
these envelopes will be enriched with Na and Al while O will be depleted due to enhanced 
extra mixing triggered by their tidal spin-up. Finally, a fraction of primordial binaries in 
globular clusters migh have been very high, up to 100% (Ivanova et al. 2005). 
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5. Concluding Remarks 

In this paper, we have elaborated upon the ideas proposed by Denissenkov & Vanden- 
Berg (2003a) about canonical extra mixing in single upper RGB stars and enhanced extra 
mixing in low-mass red giants spun-up as a result of their tidal synchronization in close 
binary systems. In order to put as many as possible observational constraints on properties 
of extra mixing, we have supplemented the old data on the Li and CN abundance changes 
along upper RGBs in field low-metallicity stars with the new data on photometry, chemical 
peculiarities and rotational periods/velocities of stars from the solar- metallicity open cluster 
M67, globular clusters 47Tuc, NGC6397, NGC6752, and RS CVn binaries. 

We have confirmed the conclusions made in the earlier paper that the secular shear 
instability driven by differential rotation of the red giants' radiative zones can be considered 
as a promising physical mechanism for both modes of extra mixing while the turbulent 
diffusion coefficient (1) derived by Maeder & Meynet (1996) can be used to model them 
appropriately, provided that: (i) unlike the Sun, low-mass MS progenitors of those red giants 
already possessed differentially rotating radiative cores; (ii) the specific angular momentum 
was conserved in each mass shell, including convective regions, inside those stars during their 
entire evolution from the MS through the RGB tip; and (Hi) the diffusion coefficient (1) is 
taken with the enhancement factor / v 20. For a binary red giant, we have additionally 
assumed that the tidal force brings to corotation only an upper part of its convective envelope. 
For the orbital and stellar parameters typical for the RS CVn binaries, this assumption results 
in spinning up of their radiative zones by a factor of > 10. 

Although we present some arguments in support of the made assumptions, it still seems 
unlikely that even our tidally enforced enhanced extra mixing with D niix ~ 10 11 cm 2 s~ x is not 
accompanied by a fast transport of angular momentum that would work toward flattening 
the f2-profile in the radiative zone, thus reducing the mixing efficiency From this point of 
view, our approach to modeling extra mixing in upper RGB stars is similar to the "Maximal 
Mixing Approach" recently used by Chaname et al. (2005). The question "how far is our 
approach from reality?" can be answered in the usual way — by comparing the theoretical 
predictions made by us with future observations. Here are some possible observational tests 
that can support or reject our models. 

First, we propose to determine 12 C/ 13 C ratios in upper RGB stars in the solar-metallicity 
open cluster NGC 6791 whose MS turn-off mass is ~ 1.1 M . If the Sun is not an exceptional 
case then, like the Sun, MS progenitors of the NGC 6791 red giants must have rotated as 
solid bodies. Therefore, we expect that canonical extra mixing in upper RGB stars in 
NGC 6791 should be inefficient. Hence, these stars should keep their post-first-dredge-up 
values of 12 C/ 13 C^ 25 unchanged instead of having the values of 12 C/ 13 C^ 10 similar to 
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those measured in the M67 evolved stars that have definitely experienced canonical extra 
mixing. 

Second, although we do believe that the depth of extra mixing is constrained by equation 
(4), when computing its rate, we still cannot choose between the constant diffusion coefficient 
(-Dmix ~ 10 9 cm 2 s _1 for canonical and D mix « 10 11 cm 2 s _1 for enhanced extra mixing) and 
the turbulent diffusion coefficient (3) that is proportional to the radiative diffusivity K. The 
first coefficient better reproduces the Na-Li data for RS CVn binaries and, unlike the second 
one, it actually results in a steeper decline of the 12 C/ 13 C ratio immediately after the bump 
luminosity as appears to be dictated by observations (e.g., compare the log 12 C/ 13 C panels 
in our Fig. 2 and in Fig. 4 from Denissenkov & VandenBerg (2003a)). Fortunately, one 
probably can discriminate between the diffusion coefficients by measuring the Li abundance 
in RGB tip and clump stars. As we have predicted, in the case of D mix oc K these evolved 
stars should have values of loge( 7 Li) > 1 comparable with the post-first-dredge-up values 
or even exceeding them. Tentatively, this theoretical prediction seems to have been already 
confirmed by Pasquini et al. (2001). However, we think more observational work should be 
done in this direction. 

Finally, we surmise that red giant components of the RS CVn binaries experience en- 
hanced extra mixing caused by their tidal spin-up. We think that the large Na abundances 
found in these stars by Morel et al. (2004) support this idea. We did not mention before 
that the same authors had also reported high Al abundances in the RS CVn binaries. It is 
interesting that a similar Na-Al correlation commonly inheres in the star-to-star abundance 
variations in globular clusters. One of its possible interpretations also employs a model of 
tidally enforced enhanced extra mixing in a binary red giant (§4.2.2). If it could be possible 
to corroborate our hypothesis of enhanced extra mixing in the RS CVn binaries observation- 
ally that would also prove that a large amount of Al can be produced at the same physical 
conditions under which Na is synthesized (this actually requires higher rates for the MgAl- 
cycle reactions responsible for the Al production). We propose to determine 12 C/ 13 C ratios 
and/or CN abundances in the RS CVn binaries to see whether they are consistent with those 
resulting from enhanced extra mixing. 
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Fig. 1. — Schematic structure of a red giant with extra mixing in its radiative zone. HBS and 
BCE stand for the "hydrogen burning shell" and the "bottom of the convective envelope" . 
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Fig. 2. — Comparison of the observational data from Gratton et al. (2000) for field metal- 
poor (—2 < [Fe/H] < —1) low-mass stars (circles) with results of our stellar evolution 
computations with rotational mixing on the upper RGB described by equation (1). It is 
assumed that the specific angular momentum is conserved in each mass shell, including 
convective regions, during the entire stellar evolution. Computations are done for M = 
0.85 M and the initial rotation parameter f £ = 0.0003 (equation (5)). Theoretical results 
are presented for the following combinations of the enhancement factor (/ v in equation (1)) 
and heavy-element mass fraction: f v = 25 for Z = 0.0005 and Z = 0.002 (solid and dot- 
dashed curves, respectively), and f v = 15 for Z = 0.0005 (dashed curve). 
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Fig. 3. — An illustration of tidal synchronization in real stars. Circles are G and K giant 
components of field binaries observed by De Medeiros et al. (2002). The factor A/n m 1.27 
takes into account the random orientation of their rotation axes. Pd is the binary orbital 
period in days. Filled circles are binaries with nearly circular orbits (eccentricities e < 0.10). 
The systems with e > 0.10 are represented by open circles. A set of theoretical solid curves 
is constructed by us assuming that Q = u, where Q and uj are the spin and orbital angular 
velocity of a red giant, for radii R/R Q = 5, 10, 20 and 40 (from left to right). The giants 
represented by filled circles are most likely to have synchronized their rotation. Dashed curve 
shows the maximum possible t> rot of a red giant that fills its Roche lobe in a binary system 
with M = 1.7 M Q , q = 0.5. 
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Fig. 4. — High precision CMD of the old solar-metallicity open cluster M67 (open circles) 
from the data of Sandquist (2004) and our 4 x 10 9 yr isochrone (solid curve). A region 
around the bump luminosity is outlined by a square. Dashed curve is the evolutionary track 
of the model star with M = 1.35 M and Z = Z Q = 0.0188. All models were computed 
with the same initial rotation parameter f £ = 0.00075. Filled circles are stars with known 
values of v sin i from Melo et al. (2001). Asterisks are red giants and clump stars for which 
the ratio 12 C/ 13 C was measured by Gilroy & Brown (1991). For the cluster reddening and 
distance modulus, we have adopted the values E(B — V) = 0.03 and (m — M) v = 9.67 that 
correspond to the observational lower limits provided by Sandquist (2004). 
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Fig. 5. — Estimates of the upper limits for the "true" mean rotational velocities of stars in 
globular clusters 47 Tuc, NGC6397 and NGC6752 by Lucatello & Gratton (2003). Filled 
square — a mean for MS turn-off stars, open square — for subgiants. Horizontal errorbars 
embrace all My values for the observed samples of stars in NGC6752 ([Fe/H] = — 1.43). 
Theoretical curve shows how t> rot changes during the evolution of our model star with M = 
0.85 M Q and Z = 0.0005 in which the specific angular momentum is conserved in each mass 
shell. 
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Fig. 6. — Observational data on (4/7r) i>sin % for M67 stars (open circles are MS stars, 
filled circles are MS turn-off, subgiant and lower RGB stars) from Melo et al. (2001) and 
theoretical dependences of v Iot on (B — V) . Solid curve — the 4 x 10 9 yr rotational isochrone 
for Z = 0.0188. Dot- short- dashed curve shows the rotational evolution of our 1.35 M model. 
Both curves are computed under the assumption of constant specific angular momentum in 
each stellar mass shell for the initial rotation parameter f £ = 0.00075. Short-dashed curve is 
the 4 x 10 9 yr isochrone with the initial velocity v rot ~ 10 kms" 1 computed assuming solid- 
body rotation inside all of the M > 1.1 M & models used to construct it. For comparison, 
dot-long-dashed curve and long-dashed curve represent the evolution of the 1.35 M model 
for the cases of differential and solid-body rotation, respectively, but for the initial velocity 
v TOt a* 40 kms -1 . 
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Fig. 7. — Profiles of some element abundances (mass fractions in panels a and b), /z-gradients 
and diffusivities (panels c and d) near the HBS in two rotating RGB models with M = 
1.35 M Q , Z = 0.0188 and the initial rotation parameter f £ = 0.00075. The models are located 
immediately above (panels b and d) and below (panels a and c) the bump luminosity. In 
panels c and d, the /i-gradients are plotted using solid curves, thermal diffusivity K — long- 
dashed curves, coefficients of the vertical turbulent diffusion D v (these are given by equation 
(1) for the enhancement factor / v = 15) — dot-dashed curves. Vertical line segments in panel 
d point to the extra mixing depths specified by parameters A log T = 0.19 (dotted line) and 
A log T = 0.22 (dashed line) as well as to the depth determined by equation (4) (solid line). 
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Fig. 8. — 12 C/ 13 C isotopic ratios in the M67 evolved stars measured by Gilroy & Brown 
(1991) (symbols with arrows and errorbars). Filled circles are RGB stars while open circles 
are clump stars. Solid curve shows a decrease of 12 C/ 13 C produced by the standard first 
dredge-up for the model star with M = 1.35 M and Z = 0.0188. Dotted curve presents 
results of our computations in which extra mixing on the upper RGB has been simulated 
using the diffusion model with the depth and rate parameters: A log T = 0.22 and D mix = 
8 x 10 8 cm 2 s _1 . Dashed curve is obtained for our rotating evolutionary model with extra 
mixing modeled by the vertical turbulent diffusion (the diffusion coefficient from equation 
(1) with the enhancement factor / v = 15). The initial rotation parameter was f £ = 0.00075. 
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Fig. 9. — Dashed curves are evolutionary tracks of single solar- metallicity stars. Their masses 
are indicated on the left from the tracks. Solid curve shows the evolution of the 1.7 M star 
placed into a binary system with q = 0.5 and a = 80 R Q . Its initial rotation is specified 
by the parameter f £ = 0.00075. After having been tidally spun up on the lower RGB, this 
binary star makes an extended bump luminosity zigzag between My m 0.7 and My m 1.5. 
Filled circles are RS CVn binaries that have nearly solar metallicities and estimated masses 
1.5 < M/M < 1.7 (Fekel et al. 2002; Morel et al. 2004; Fekel & Henry 2005). Each of them 
has almost equal rotational and orbital periods that correspond to a ~ 50 R® or a ~ 80 R Q . 
Open circles are binary red giants with circularized orbits from Fig. 3. Asterisk is the star 
HD 21018 that lies above the dashed curve in Fig. 3. 
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Fig. 10. — Dot-long- dashed and dot-short- dashed curves are the ZAMS and bump luminosity 
f2-profiles in the single model star with M = 1.7 Mq, Z = 0.0188 and initial rotation pa- 
rameter f £ = 0.00075 rotating differentially (with the specific angular momentum conserved 
in each mass shell). Long-dashed and short-dashed curves are the ZAMS and bump lumi- 
nosity fi-profiles in the same model but for the case of its evolution with uniform rotation 
starting on the ZAMS with v rot ~ 100 kms -1 . Solid line segments show the orbital angular 
velocities lo 1 = 0.2834 and uj 2 = 0.1400 (in units of 10 _5 rads _1 ) for binary systems with a 
1.7 M e primary component, the same mass ratio q = 0.5 but with different semi-major axes 
a = 50 Rq and a = 80 R & , respectively. Vertical dotted line marks location of the BCE in 
the bump luminosity models. 
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Fig. 11. — Bump luminosity and a later f2-profiles in differentially rotating star with M = 
1.7 M , Z = 0.0188 and initial rotation parameter f £ = 0.00075. Dot-dashed curve — in a 
single star. Solid curve — in a binary star with q = 0.5 and a = 80 R Q . It is assumed that 
only an upper part of convective envelope is spun up by the tidal force (for details, see text). 
Short-dashed curve — a later profile evolved from the solid one. Long-dashed curve — in the 
same binary star but under the assumption that rotation of the whole convective envelope 
has been tidally synchronized. Dotted curve — a binary star with q = 0.5 and a = 50 R Q . 
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Fig. 12. — Theoretical dependences of the surface rotational velocity t> ro t on My for the 
solar-metallicity 1.7 M primary star in binaries with the same q = 0.5 but different semi- 
major axes: a = 50 Rq (dashed and dotted curves) and a = 80 R & (solid and dot-dashed 
curves). These are compared with locations of RS CVn binaries (circles are data from Fekel 
et al. 2002; Morel et al. 2004; Fekel & Henry 2005). Solid and dashed curves start with a 
subgiant model having R = 3 R Q and Q = 0, while dot-dashed and dotted curves initially (at 
R = 3R Q ) have f2_ 5 = 2.5, which extrapolates to v rot ~ 100 kms -1 back on the ZAMS. The 
evolution of v rot = Rfl and a due to the tidal interaction is followed by solving equations 
(7-8) as described in text. 
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Fig. 13. — Tidal evolution of the ratio of the spin and orbital angular velocities for a solar- 
metallicity 1.7 M primary star in a binary system with the same q = 0.5 but different 
semi-major axes: a = 50 R & (dashed curve) and a = 80 R Q (solid curve). We define the 
corotation radius R COT as the radius at which Q/u = 0.9. During the subsequent evolution 
of the primary star it is assumed that the part of its convective envelope above R cor rotates 
as a solid body with Q = u, while the part at -Rrce < r < Rcor still rotates differentially. 
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Fig. 14. — The single and binary evolutionary models of a solar-metallicity 1.7 M star with 
the corotation radii R cor /R Q = 2, 3, 4, 5, 6, 7, and 8 (from bottom to top). The increase of 
R cor causes the radiative zone to spin up (panel a), which makes effects of rotation on the 
model's internal structure to be more pronounced, resulting in a noticeable extension of the 
bump luminosity zigzag (panel b; in order not to overlap each other, the tracks are shifted 
along the ordinate). 
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Fig. 15. — Differential luminosity functions constructed using the evolutionary track of the 
single solar-metallicity 1.7 M & star (thick curve) and the track of the same star placed into 
the binary system with q = 0.5 and a = 80 R Q (thin curve). The time spent by the single 
star between My = 2 and My = is 109 million years, while for the binary star it takes 
146 million years to make an extended zigzag on the CMD when My increases from 0.71 to 
1.52 and then decreases back to 0.71. For the binary star, the bump luminosity is shifted 
towards the lower RGB by AM V m 0.8. 
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Fig. 16. — Profiles of some element abundances (panels a, c) and diffusivities (panels b, d) 
in the single solar- metallicity 1.7 M Q model star (panels a, b) and in the same model placed 
into the binary system with q = 0.5 and a = 80 R Q (panels c, d). Both models are beyond 
the bump luminosity. In panels b and d, dot-dashed curves show the turbulent diffusion 
coefficient (1) for / v = 1, dashed curves — thermal diffusivity, solid curves — /z-gradients, 
and dotted curves — constant diffusion coefficients for canonical (D m[x = 10 9 cm 2 s _1 , panel 
b) and enhanced extra mixing (D m[x = 10 11 cm 2 s _1 , panel d). The mixing depth is always 
constrained by equation (4). 
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Fig. 17. — Peculiar combinations of the Li and Na abundances in RS CVn binaries (circles), 
as reported by Morel et al. (2004), and their theoretical reproduction with our model of 
enhanced extra mixing in the solar-metallicity 1.7 M & binary red giant. Dashed curve is 
obtained using the diffusion coefficient (1) with f v = 20, solid curve — for the constant co- 
efficient -D m ix = 10 11 cm 2 s -1 with the mixing depth from equation (4). Dotted line segments 
delimit regions of the bump luminosity zigzag extended by fast rotation. 
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Fig. 18. — The O-Na anticorrelation in red giants in the globular clusters M3 (circles) and 
M13 (triangles) from a recent update by Sneden et al. (2004). Filled symbols — stars with 
lower surface gravities (higher luminosities). Solid curve shows a simple mixture of a (1 — in- 
fraction of material with the abundances [0/Fe] ini t = 0.3, [Na/Fe] init = —0.2 and a x-fraction 
of material with [0/Fe] acc = —1.0, [Na/Fe] acc = 0.7 (thick dots on the curve mark values of 
x from to 1 with the increment 0.1). Dashed curve is obtained with our model of tidally 
enforced enhanced extra mixing in a red giant with M = 1.0 M©, Y = 0.24, Z = 0.0005 
and f £ = 0.0003 in a binary system with q = 0.3 and a = 50 R®. It is computed using the 
diffusion coefficient (1) with / v = 20. Dotted curve presents a test case in which Q in the 
radiative zone has been increased by the factor of 25 as compared to the single star case. 
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Fig. 19. — The same as in Fig. 18 but for the Na-Al correlation. For Al, the initial and 
most extreme abundances are [Al/Fe]; n i t = —0.2, [Al/Fe] acc = 1.2. Following Denissenkov & 
Tout (2000), we have additionally assumed that [ 25 Mg/Fe]i n i t = 1.2 and the rate of reaction 
26 Al g (p,7) 27 Si is a factor of 10 3 as fast as its rate from Angulo et al. (1999). 



